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We propose an optical scheme to generate a superposition ol coherent states with enhanced size 
adopting an interlerometric setting at the single-photon level currently available in the laboratory. 
Our scheme employs a nondegenerate optical parametric amplifier together with two beam splitters 
so that the detection of single photons at the output conditionally implements the desired superpo- 
sition of second-order photonic operations. We analyze our proposed scheme by considering realistic 
on-off photodetectors with nonideal efficiency in heralding the success of conditional events. A 
high-quality performance of our scheme is demonstrated in view of various criteria such as quantum 
fidelity, mean output energy, and measure of quantum interference. 

PACS numbers: 42.50.Dv, 42.50.Xa 



I. INTRODUCTION 

Since Schrodinger envisioned a Gedankenexperiment to 
illustrate some peculiar feature of quantum mechanics 
when applied to macroscopic objects [TJ, the so-called 
Schrodinger cat state has long been of great interest 
from a fundamental point of view. Furthermore, the 
Schrodinger cat state has also been identified as a use- 
ful resource for practical applications, e.g., in the field 
of quantum information science over the past decades. 
There have been a great deal of theoretical and ex- 
perimental attempts to generate a Schrodinger-cat-type 
state in various physical systems including atomic system 
[3] , solid-state system [3] , and optical system [SHU] • 
Among these, an optical-cat-like state, especially a su- 
perposition of coherent states (SCS) [T31 H3] , turns out 
to be useful for quantum computation [T5Tf2"0] . quantum 
teleportation , quantum-enhanced metrology [24]— 

I29j . and fundamental test of quantum mechanics |30l - 
139] . Along this line of practical significance, a number 
of proposals have been made to generate an optical SCS 
[TBI I40H5T] , which subsequently led to considerable ex- 
perimental progresses [6HI2], 

The original proposals for generating an optical SCS 
in a deterministic way were based on Kerr nonlinear 
medium |13l I40j . However, their realizations have not 
yet been made since the high nonlinearity required for 
such a deterministic scheme is not achievable with the 
current technology. Alternatively, most of the SCSs cre- 
ated so far in laboratory [BHU] are based on a nonde- 
terministic scheme that generates a target state only un- 
der prespecified conditions. For instance, the photon- 
subtraction method |42j can implement a SCS relying 
on the fact that even (odd)-photon-subtracted squeezed 
vacuum states well approximate even (odd) SCSs with 
small amplitudes. Marek et al. also investigated the 
effect of photon addition as well as subtraction on the 
initial squeezed vacuum state to obtain an output state 



close to a squeezed SCS [50]. On the other hand, the ex- 
perimental scheme of Ref. [5] produced a squeezed SCS 
by conditioning an input Fock state on the outcome of 
homodyne detection. That is, if a Fock state |n) is in- 
jected into a beam splitter and conditioned on a certain 
predetermined value of the quadrature amplitude at one 
output mode, the other output state converges well to a 
3-dB-squeezed SCS as n becomes large. 

In this paper we propose an experimental scheme to 
generate a squeezed SCSs with enhanced size employ- 
ing a coherent superposition of photonic operations. Our 
proposal makes use of an interferometic setting that com- 
bines a nondegenerate parametric amplifier (NDPA) to- 
gether with two beam splitters. The detection of sin- 
gle photons at two output modes conditionally imple- 
ments a coherent superposition of second-order opera- 
tions, namely, aa) and a) 2 , on an arbitrary input state. 
We show that this scheme, when applied to a single- 
photon input, can successfully generate an output state 
very close to an odd-parity squeezed SCS which includes 
a triple-photon component or higher. To illustrate the 
feasibility of our scheme, we consider the use of nonideal 
on-off photodetectors for heralding the conditional events 
and demonstrate a high-quality performance in terms of 
various criteria such as quantum fidelity, mean output en- 
ergy, and the measure of quantum intereference [52j . Our 
proposed scheme thus seems very feasible within the cur- 
rently available experimental techniques [53] . In particu- 
lar, we note that a single-photon interferometric scheme 
similar to ours was recently proposed [54] and also exper- 
imentally realized [55] to verify the bosonic commutation 



relation [a, fit] = 1. 

This paper is organized as follows. In Sec. [n] we briefly 
introduce a class of squeezed SCSs targeted in our work. 



In Sec. Ill we propose an experimental scheme to imple- 
ment a superposition operation to our end. With a brief 
account of the working principle under our scheme, its 
rigorous description is also provided within the Wigner- 
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function formalism. In Sec. IV we evaluate the per- states turn out to be 



formance of our scheme in terms of various criteria such 
as quantum fidelity, mean output energy, and measure 
of quantum interference. In Sec. [Vj our analysis is fur- 
ther extended to investigate the effect of on-off photo de- 
tectors with nonideal efficiency. Finally, our results are 
summarized in Sec. I VII 



II. SQUEEZED SCS 

Although there exist some generalized versions of SCSs 
[131 Ej, we begin with two typical SCSs, namely, even 
and odd SCSs, 



|SCS±(a)) =M ± (\a)±\-a)), 



(1) 



where Af± = 1 / 'v '2 ± 2e" 2 " 2 are the normalization fac- 
tors and |±a) is a coherent state of amplitude ±a. Here 
a is assumed to be real for simplicity but without loss of 
generality. The SCS with the positive (negative) sign in 
Eq. ([I]) is called an even (odd) SCS because it contains 
only even (odd)-numbered Fock states for any a. The 
size of SCS may be quantified by the magnitude of the 
amplitude |a|, which also characterizes the distance be- 
tween the two component states \a) and \—a) in phase 
space. 

Suppose now that we apply a squeezing operation 



S(s) 



exp [f 



)] (s: real) on the SCS in Eq. 



11). Note that S (s) squeezes a quantum state along the 
real (imaginary) axis in phase space for a positive (neg- 
ative) s. When applied to |SCS±) in Eq. |l]), S (s) does 
not change the parity of |SCS±) since S (s) is invariant 



under the parity operator 
the state 



-l) a °. Hereafter, we label 



3 SCS±(a)) =S(«)|SCS±(or)) 



(2) 



as a squeezed (even or odd) SCS (sSCS). In view of Eq. 
([2]), let us now define a class of states, 



|^ n )=,S(lnV2)|SCS ( _ ) „(V^)>, 



(3) 



which is 3-dB-squeezed SCS and has the parity (—1)" 
and amplitude a = i/n for a non-negative integer n. 

Ourjoumtsev et al. proposed a novel conditional 
scheme to create a similar Shrodinger-cat-type state, 
which generated a state close to \ip2) in Eq. ([3]). Their 
idea was to utilize homodyne-conditioning technique with 
ri-photon state \n) as an input [5]: The state \n) is in- 
jected into a 50:50 beamsplitter with the other input in 
a vacuum state. The output state is kept conditioned on 
a particular outcome of homodyne detection at the other 
mode, which produces a class of states \<f> n ). Along with 
the trivial states \<f)o) = |0) and \4>i) = other lowest 




(4) 
(5) 
(6) 



It is worth noting that \4> n ) has the same parity as the 
Fock state \n): If n is even (odd), \cj> n ) has only even 
(odd) number of photons just like |SCS±) and |sSCS±). 

Interestingly, \cj) n ) is very close to \ip n ) in view of quan- 
tum fidelity: F = \ (ip n |0„)| 2 > 0.97 for all n > 0, and 
F — > 1 as n -> 00. Ourjoumtsev et al. produced a 
state targeting |^ 2 ) [SS] with a two-photon state |2) as 
an input. The target state Ifc) is very close to |^ 2 ) with 
the fidelity F = 97.2% and can be made even closer to 
|sSCS±) with F = 99.0% by increasing the squeezing 
level from 3 to 3.5 dB. In the next sections, we propose 
an experimental scheme to generate a larger-size sSCS, 
i.e., \(f>3,) and beyond. 



III. SCHEME TO GENERATE SCS 

In this section, we present an experimental scheme to 
implement a coherent superposition of photonic opera- 
tions for generating a sSCS. For instance, to generate 
the target state ^3), which is a superposition of the 
number states |1) and |3), one has to apply a super- 
position of identity operation / and two-photon addi- 
tion a^ 2 to an input single-photon state |1). That is, 
(tl + ra^ 2 )|l) ~ ci|l) + C2 [ 3) . Our scheme makes use 
of a NDPA together with two beam splitters, and the 
output state is kept on the condition of single-photon de- 
tections at two output modes of a beam splitter (Fig. 2). 
This implements the superposition operation of the form 
taa^ + ra) 2 , as will be shown below. Note that the effect 
of the operation aa^ is essentially the same as that of the 
identity operation / for the case of Fock-state input. 

To begin with, we describe our experimental scheme 
under the small-parameter approximations |54|, I55[ I57j , 
which characterizes the output state when using an 
NDPA with small squeezing. In Sec. Ill B, we will more 
rigorously describe our proposed scheme using a Wigner- 
function formalism in phase space. 



A. Brief sketch of the scheme 

The first stage of our scheme is to prepare a single- 
photon state |1). This can be done by pumping an NDPA 
with both inputs of signal and idler modes in vacuum 
states and then by projecting the output idler to a single- 
photon state (Fig. [l]). More precisely, the interaction 
between the signal mode a and the idler mode d within 
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FIG. 1: (Color online) Experimental scheme for a prelimi- 
nary single-photon generation. The pump beam creates twin 
beams as it passes through a NDPA. One output mode is 
directed to the detector PD,j and the other to a mirror. A 
detection event at PDd heralds the generation of an approxi- 
mate single-photon state in mode a. 



the NDPA is described by a two-mode squeezing operator 
S ad (s)=exp[s(ad-aU% (7) 

With the two inputs in vacuum states, the output be- 
comes 



m ad =s ad (s)\o) a \o) d 

« (l - B tf$) \0) a \0) d 
= |0)JO) rf - S |l)Jl 



d • 



(8) 
(9) 



where the condition s <C 1 (small squeezing) is used. 
Now, the detection of a single photon at mode d gives 
the desired single-photon state in mode a: 



l*> Q 



<i| d (l*U 



ID 



(10) 



In Fig. [T] the signal mode a is reflected by a mirror 
to be injected again into the same NDPA (see Fig. 2). 
By doing this, we can reduce the number of required 
NDPAs, and furthermore, as the input signal has been 
produced within the same NDPA, the mode-matching 
condition, which is critical for quantum interference, may 
be better satisfied. Thus, in the next stage, we reuse 
the pump beam during its coherence time for the NDPA 
process. Under the small-squeezing condition, the output 
now reads 

l*) ac =^oc(*)|*) a |0) c (11) 

« (l-^ct + i s 2 at 2 ct 2 )|l)JO) c (12) 
= |1)J0) c -^ S |2)J1) c + V3 S 2 |3)J2) c . (13) 

Note that other irrelevant terms are not included in Eq. 
( 12 1 and that only the last two terms in Eq. (131 will sur- 



vive the subsequent detection process as described below. 

In the third stage, we tap a very small amount of the 
beam in mode a using a beamsplitter BS a {, whose trans- 
mittance is very close to 1. This can be described by a 




output 



FIG. 2: (Color online) A whole process to generate a de- 
sired state |(^3} and beyond. After generating the approxi- 
mate single-photon state, the pump beam is redirected by a 
mirror into the NDPA and creates another twin beams. The 
coincident detection of single photons at PDf, and PD C her- 
alds a successful implementation of a superposition operation 
coherently combining aa^ and a' 2 . See the main text. 



beamsplitter operator B ab {t) acting on the state | v E') ac , 
with mode b in a vacuum state, as 



b ' 



where 



Babit) = cxp 



(a6 f - a r bj 



(14) 



(15) 



with its transmissivity (reflectivity) given by t = cos | 
(r = sin |). Under the condition r± -C 1 (ti ~ 1) and 



0{r\) ~ O(s), Eq. (14 1 may be approximated as 



l*>«6c« (l + ^ t )l*) M |0)6 (16) 

~ -2sn \l) a \l) b |l) c + V3s 2 \3) a |0) b |2) c .(17) 

Finally, after the state \^} abc passes through another 
beam-splitter BS& C (transmissivity t), we obtain the out- 
put state 

m'abc = B bc (t)\*) abc (18) 

«2 S r 1 (r 2 -t 2 )|l)Jl)Jl) c 

+V3sHr |3)J1)J1) C . (19) 



Note that the first term in ( 19 ) is due to the first term in 
( IT): The factor i 2 (r 2 ) indicates that two single photons 
are simultaneously transmitted (reflected) via the BSf, c , 
and the different sign between r 2 and t 2 arises due to the 
phase change induced by two reflections at the beam- 
splitter. The last term in ( 19 ) implies that two photons 



of mode c in (17) are split into individual single photons 
in mode b and c. The whole process here is a quantum- 
mechanical superposition of three different paths, which 
is illustrated by the diagrams in Fig. [3j 
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[2 S ;;{r-r)|l) a+ ^r|3)jl)Jl) £ 

FIG. 3: (Color online) The first two and the third diagrams 
show a single- and triple-photon generation, respectively. Yel- 
low arrows represent the single-photon generated in the first 
stage of our scheme in Fig. 2, whereas red arrows represent a 
single-photon in the first and the second diagram, or a two- 
photon in the third diagram, generated in each mode b and c. 
The parameter s denotes the degree of squeezing within the 
NDPA, n the reflectivity of the beamsplitter BS a i,, and t (r) 
the transmissivity (reflectivity) of the beamsplitter BS(, C - In 
the first two diagrams, s indicates that twin single-photons 
are generated and n indicates that upper one of twin beams 
is reflected by BS a (,. Two t's (r's) in the first (second) dia- 
gram means that two beams all pass through (are reflected at) 
BStc without reflection (transmission). In the third diagram, 
s 2 indicates that twin two-photons are generated and t and r 
means that one of the two photons in the lower beam passes 
through, and the other is reflected, at BSj, c . 



Now, if both detectors PD D and PD C simultaneously 
register single photons, the conditional output is given 
by 



B. Rigorous formulation using Wigner- function 
approach 

In Sec. |III A| we have briefly described the work- 
ing principle of generating the target sSCS under small- 
parameter approximations. We now present an exact for- 
mulation of our experimental scheme in Fig. 2 using the 
Wigner-function formalism in phase space. We here give 
a concise description of the Wigner-function method with 
essential steps and all detailed derivations are shown in 
Appendix B. 

Let us start with the twin-beam state created within 
the NDPA, which is formally described by the two-mode 
squeezing of the input vacuum states |0) o \0) d . In phase- 
space representation, a two-mode squeezing has the ef- 
fect of pscudorotation on the arguments of two quasi- 
distribution functions. That is, the input two-mode 
Wigner function W vac ((3 a )W vac ((3 d ) transforms into 

W ad (Pa, (3d) = Wvac (Pa COsh S + (3*. sinh s) 

x Wvac {(3 d cosh s + (3* sinh s) 
= (2/n) 2 exp [ - 2 cosh 2s( \/3 a \ 2 + \/3 d \ 2 ) 

+2 sinh 2«()9 fl, ( 22 ) 

where W vac {(3) = (2/tt) exp(-2|/3| 2 ) is the Wigner func- 
tion of a vacuum state |0), with (3 — (3 r + ifa. 

If the detector PD^ in Fig. 2 registers any number of 
photons, the conditional output state is given by p a ~ 
Ti'd [(/ — |0)(0|)d Pad], whose Wigner function becomes 

W a (J3 a ) = (l/P d ) J d 2 (3 d W ad (/3 a ,Pd) 

x [l-nW vac (p d )}, (23) 



l*U~<l|6<l|c(l*>a6c) 

« 2sn(r 2 - t 2 ) \l) a + V3s 2 tr |3) a . (20) 

Thus, by adjusting the parameters s, ti(ri), and t[r) in 
the experimental scheme of Fig. 2, we can obtain 1 03^ in 
Eq. (§. 

More generally, if an arbitrary input state |^) in is in- 
jected into the scheme in Fig. 2, one can similarly confirm 
that the output state is given by 

Wout = (lib (lie B bc (t)B ab (h)S ac |*) in |0) 6 |0) c 

w [C x aa) + C^ 2 ) |*) in (21) 



where C\ = sr\[r 2 — t 2 ) and C% — -^s 2 tr. There- 
fore, our scheme in Fig. 2 implements a superposition 
of second-order operations aa' and a) 2 on an arbitrary 
input. In principle, one can generate \4>§) , \4>b) , \<f>r) , ••■ 
or \4>2} , |<^>4) , |06) !••• by consecutively applying this su- 
perposition operation to an input state |1) or |0), respec- 
tively. 



with its normalization 

P d = Tr ad [I a ®(I- \0){0\) d pad] 

d 2 (3 a ld 2 p d W ad (P a , (3 d ) [1 - nW mc {(3 d )] (24) 



Note that P d represents the conditional probability for 
detection events and its detailed funtional form is given 
in Appendix: B. Equation (23) is an exact version of Eq. 



(101 



The conditional state in (23 1 is again injected to the 



NDPA with other mode c in a vacuum state. We thus 
obtain an output state from the NDPA, whose Wigner 
function is given by 

W ac (p a , (3 C ) = W a (fi a COsh S + (3* c sinh s) 

x Wvac ((3 C cosh s + (3* a sinh s) . (25) 

The third step, i.e., mixing the mode a with another 
mode b (vacuum state) at a beamsplitter B a b(t 2 ~ 1 or 
r\ <C 1) is described by 

W abc (p a , (3 b , (3 C ) = W ac (hp a + rift, (3c) 

x W vac (h/3 b - n0 a ) . (26) 
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(a) 



(b) 



(c) 



FIG. 4: (Color online) Wigner functions of (a) \4>3), (b) a state 
generated via the scheme in Fig. 2, and (c) sSCS [Eq.|2|]. 
The state in (b) is produced with the parameters s — 0.04 
and t{ = 0.999 (r? = 0.001). The sSCS in (c) is 2.9-dB 
squeezed (s' = 0.33) with a — 1.7. The fidelity of the output 
state in (b) is 99% with respect to the sSCS. 



Another mixing process at the beamsplitter B bc gives 

W' abc {P a ,P b , Pc) = W abc (p a , tp b + rp c , tp c - r/3 b ) . (27) 

Finally, the success probability for the coincident detec- 
tion of photons at PDb and PD C is given, similar to Eq. 
PI, by 



P bc = J d 2 p a J d 2 p b J d 2 p c W abc (p a ,f3 b ,f3 c ) 

x [1 - ttW v&c (&)] [1 - 7rW vac {p c )\ . (28) 

In addition, the Wigner function of the conditional out- 
put is given by 

Wo ut (P a ) = (l/P bc ) Jd 2 p b J d 2 P c W' abc {Pa, Pb, Pc) 

x [1 - nW mc (P b )} [1 - irW mc (P c )] . (29) 

One of the typical examples generated under our 
scheme is shown in Fig. [3J together with the target state 
\4> 3 ) and sSCS [Eq.pJ]. The figures show that the output 
state is slightly closer to sSCS than to \4>3). We investi- 
gate the quantum fidelity with respect to sSCS as 



(30) 



F = n d z p W out {p) W sSGS (P) ■ 



where VV" sSC s (P) is the Wigner function of sSCS [Eq.Q] 
with squeezing parameter s': 

Wscs(P) =N 2 [W mc (l3r - a, Pi) + W V3JC (P r + a, Pi) 

-2W vac (/3)cos(4aA)]. (31) 
W sS cs(P) = W SCS (Pcoshs' + p* sinhs') . (32) 

In Fig. 4, the quantum fidelity of the output state is 99% 
with respect to the sSCS. 

IV. PERFORMANCE OF THE SCHEME 

In this section, we investigate the overall performance 
of our scheme in terms of various criteria. In particular, 
we look into success probability, quantum fidelity, mean 
photon number, and a measure of quantum interference 
to properly address the quality of an output state under 
our scheme. 
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FIG. 5: (Color online) (a) success probability and (b) quan- 
tum fidelity of an output state as functions of the squeezing 
parameter s. The three curves in each panel represent the 
casus of r'l = 0.001 (blue solid), 0.01 (red dashed), and 0.1 
(green dot-dashed), respectively. The success probability P 
is the joint detection probability of all three detectors in the 
scheme of Fig. [2] The fidelity F quantifies the closeness 
between the generated sta te and the target sS CS, which is 
optimized over t (r) in (271, s' in (32 I and a in (31 1. 



We first investigate the success probability P and the 
fidelity F in Eq. (30) in more detail. The overall suc- 
cess probability P is given by P — PdP bc , which is the 
joint probability of generating a single-photon state in 
the first stage and detecting photons simultaneously at 
PDf, and PD C in the final stage. Figure [5] shows the 
success probability and the fidelity as functions of the 
squeezing parameter s. In each panel, three curves refer 
to the cases with different reflectances of the beams plit- 
ter BS ab in Fig. § i.e., rf = 1 - 1\ = 0.001 (blue solid), 
0.01 (red dashed), and 0.1 (green dot-dashed). Note that 
the fidelity is optimized over s', a, and t (or r) at fixed 
s and ti and that the success probability is evaluated 
corresponding to those parameters. 

Some general trends can be seen from Figs. [5^a) 
and[5]jb). First, the success probability increases with 
the squeezing parameter s whereas the quantum fidelity 
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decreases with s. Second, similarly, increasing the re- 
flectance r\ enhances the success probability but lowers 
the fidelity as a trade-off. A higher squeezing produces 
more twin-beam photons within the NDPA, hence more 
chances of detecting photons at the detectors PD& and 
PD C . Similar reasoning can be given to the trend with 
the reflectance of BS a f,. On the other hand, the degra- 
dation of quantum fidelity with s and r% is somewhat 
expected as our scheme is envisioned under the condi- 
tion that the squeezing s is small and that r\ should also 
be small in the same order of s, as discussed in Sec. |III A| 

Even though the quantum fidelity can be very high, it 
would be practically not useful if the conditional proba- 
bility is too low. However, in our scheme, if the fidelity 
above 90% is acceptable, the success probability goes up 
to 1CT 6 - 1CT 5 for r\ = 0.001 or 0.01. Thus, if the rep- 
etition rate of a pulsed beam is assumed to be around 1 
MHz, the generation rate becomes 1-10 Hz. Therefore, 
with moderate values of s and n, the present scheme 
seems quite feasible for an experimental realization of 
high-fidelity sSCS. 

Next, as a size criterion, we investigate the mean pho- 
ton number (n) of the generated state. As can be seen 
from Fig. p[a), a larger squeezing s and a smaller re- 
flectance rfgive a larger size of the generated state, which 
may be expected from the relative ratio of the component 
states in Eq. (20). Note that the mean photon number 
(n) can go over 3, which is relatively large compared to 
2.75 of the recently generated state [12]. However, as 
can be seen from Fig. [6] (b) , in order to obtain a larger 
size state, the fidelity should be sacrificed to some de- 
gree. Nevertheless, with the size being equal to 2.75, the 
fidelity 83% under our scheme is still higher than 59% in 
Ref. [12]. 

Finally, we consider another kind of size criterion to 
assess the quality of quantum superposition, which was 
very recently proposed by Lee and Jeong [52j . Their mea- 
sure provides a well-defined criterion for superposition 
states by quantifying the potential of quantum interfer- 
ence of a state, which can be visualized in phase space. 
For this reason, hereafter, we call this measure the "mea- 
sure of quantum interference" (MQI) and define it for a 
single mode case as 



(a) 



MQI - | Jd 2 (3W((3) 



d 2 



d/3df3* 



1 



W(f3) , (33) 



in terms of Wigner function W(f3). A phase-space for- 
malism using a quasi-probability-distribution addresses a 
general quantum state, so our measure can be used not 
only for a pure superposition state but also for a mixed 
state. Further, MQI quantitatively manifests unique 
properties of superposition states. It gives maximum val- 
ues for pure superposition state and nearly zero values for 
fully mixed states. It can be proved that the MQI of a 
state is bounded by its mean particle number (n) and 
reaches this maximum value only if it is a pure super- 
position state. Consequently, MQI provides useful and 
sensible information on a superposition state in regard 
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FIG. 6: (Color online) (a) Mean photon number (n) as a 
function of the degree of squeezing s for the cases of r\ = 0.001 
(blue solid), 0.01 (red dashed), and 0.1 (green dot-dashed), 
respectively. Here the parameters used for optimizing the 
fidelity in Fig. [5] are fixed as t = 0.21, s' = 0.57, and a = 
2.4. Note that (h) increases with s but decreases with the 
reflectance r\. Since a larger squeezing degrades the fidelity 
see Fig. [BJb)], the fidelity as a function of (n) is also plotted 
in (b). With a fixed target state, there is a critical value of s 
which gives the maximum value for the fidelity. Thereafter, 
(h) increases, but the fidelity decreases, with s. 



to the size of its coherence. For the present scheme, we 
attempt to assess the quality of generated states by eval- 
uating the ratio of MQI to the mean photon number (n) . 
This quantity may provide an insight to the relative size 
of "superpositionness" of a quantum state. It would give 
a value close to unity if the probed state is close to a 
pure superposition state and decrease as the state loses 
its coherence and hence becomes more mixed. 

Figure [7] shows the original and the relative values of 
MQI for the states considered in Fig. [6j As can be seen 
from these plots, the value of MQI increases with the 
squeezing s for some cases. However, this is mainly due 
to the increased size of the state, not due to the increased 
coherence. This is confirmed in the plot of relative MQI 
see Fig. [7JM], where one can see that the relative values 
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FIG. 7: (Color online) (a) MQI and (b) relative MQI (MQI 
divided by (n)) as functions of the degree of squeezing s for 
the case of r\ — 0.001 (blue solid), 0.01 (red dashed), and 0.1 
(green dot-dashed), respectively. 

of MQI monotonically decrease with s. Note that the 
qualitative behavior of the relative MQI in this plot is 
very similar to the fidelity in Fig. [5] In this sense, the 
relative MQI seems to provide another useful criterion for 
the quality of a superposition state as well as its fidelity. 
Moreover, MQI has an advantage over the fidelity such 
that it can assess the quality of a generated state without 
reference to a target state. 

V. PRACTICAL APPLICATION 

In this section, we further investigate the experimental 
feasibility of our scheme by considering the inefficiency 
of devices and then comparing its performance with the 
one in Ref. [T5]. First, we take into consideration the in- 
efficiency of detectors, which can be modeled by placing 
an imaginary beamsplitter in front of each detector and 
varying its transmittance, say rj. For a perfect (totally 
inefficient) detector, 77 = 1 (0). For a low success probabil- 
ity, dark counts may also play another detrimental role. 
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FIG. 8: (Color online) (a) The success probability P and 
(b) fidelity F as functions of detector efficiency r\ at squeez- 
ing s = 0.16. The three curves denote the cases of r? = 0.001 
(blue solid), 0.01 (red dashed), and 0.1 (green dot-dashed) re- 
spectively as before. Notice that as r\ decreases, P and F also 
decrease but that, nevertheless, the fidelity is relatively robust 
against the inefficiency compared with the case of probability. 



However, it is reported that a coincidence scheme — e.g., 
which records only the synchronized outcomes of laser 
pulse signal and a click of a detector — significantly re- 
duces dark count events [20] • There has also been a re- 
port of recent progress on on-off type detectors with high 
efficiency and negligible dark count rate [55] . 

We again calculate the success probability P and the 
fidelity F, but now as functions of efficiency r\ for all the 
detectors used in our scheme. We plot P and F for the 
case of moderate squeezing in Fig. [8] At s = 0.16, the 
fidelity for r\ = 0.001 (rf = 0.01) is approximately 90% 
(89%). Generally, P and F decrease as the detector be- 
comes less efficient, i.e., with r\ decreasing. However, de- 
spite the rapid drop of the probability, the fidelity is very 
robust against the detector efficiency and still remains 
high. Moreover, this robustness of the fidelity can also 
be demonstrated together with the mean photon number. 
With the experimental report of Ref. |12j in mind, we 
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FIG. 9: (Color online) Fidelity F and mean photon number 
(n) for the cases of r? = 0.001 (blue solid), 0.01 (red dashed) 
at r) = 1, 0.9, 0.1. Refer to Fig. |6|b) for the starting val- 
ues of F and (n) at 7] = 1. As expected, F decreases with r\ 
decreasing, while (n) slightly increases. However, the degra- 
dation of the fidelity is not very significant in both of the cases 
r? = 0.001 and 0.01. 



investigate the effect of detector efficiency i] for the case 
of (n) = 2.75 in [H] and the fidelity F = 83% (74%) at 
r\ = 0.001 (0.01) in Fig. ^b). Figure [9] shows both the 
fidelity and the mean photon number as rj decreases from 
1 to 0.1 by a step of 0.1. As can be seen from the plot, F 
is degraded with a decreasing 77 while (n) rather increases 
in the opposite way. However, the degradation of F is not 
very significant in both the cases of r\ = 0.001 and 0.01 
and the fidelities are still high compared to the one in [T3] . 
Even when 77 is lowered to as small as 0.1, the fidelities are 
still higher than 59% reported in [12], with mean value 
slightly higher than (n) = 2.75. Further, if we only fo- 
cus on maximizing the mean photon number, our scheme 
can generate a sSCS as large as (n) — 3.24 (3.18) with 
predetermined squeezing s = 0.2 and fidelity F = 59% 
for r\ = 0.001 (0.01) even though the detectors have the 
detector efficiency as low as 77 = 0.1. Finally, we also 



present in Fig. 10 the moduli of the density-matrix en- 
tries of the generated state with the squeezing s = 0.16 



0.001. Figure 10 shows that the single- and 



and r\ 

triple-photon components dominate the other ones and 
that even-number components do not appear even when 
the detector efficiency is very low. This means that our 
scheme is reasonably robust against the detector ineffi- 
ciency and can generate definite odd-parity states even in 
the presence of imperfect detectors. These aspects thus 
illustrate the usefulness of our scheme for generating a 
viable sSCS if implemented in a real experimental setup. 



VI. CONCLUDING REMARKS 

We have studied an experimental scheme to generate 
and grow a large-size Schrodinger-cat-type state using an 
interferometric setting that implements a coherent su- 
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FIG. 10: (Color online) The moduli of the entries of the den- 
sity matrix of the output state which is generated with the 
parameters s = 0.16, r\ = 0.001, and the detector efficiency 
(a) 17 = 1, (b) 77 = 0.5, (c) 77 = 0.1. 



perposition of the operations aa^ and a* 2 on an arbi- 
trary input state. The proposed scheme makes use of an 
NDPA together with two beam splitters, and the detec- 
tion of single photons at two output modes conditionally 
realizes the superposition operation of our interest. We 
investigated the performance of our scheme by employ- 
ing various criteria such as the success probability, the 
quantum fidelity, the mean output energy, and the mea- 
sure of quantum interference. We have shown that our 
scheme can generate an output state of high quality in 
view of all these criteria under the use of nonideal on- 
off detectors for heralding conditional events. Therefore, 
together with the experimental achievements reported re- 
cently [12] [55] , it seems that our proposed technique can 
be useful for generating a large-size nonclassical super- 
position state within existing technology. 



VII. APPENDIX: A 

In this section, we give detailed analytical forms to the 
classes of states defined in the main text. For example, 
the state |sSCS±) is represented by its wave- function as 



(x IsSCS 



A4 



- ; 



± exp 



2 s 



exp 



V2, 



ae 



V2c 



2e" 2s 
\ 2- 



2e- 



(34) 



valid up to an overall phase [ 
state of the quadrature operator x 
detailed form of \x) is given by [55] 



Here \x) is the eigen- 



(a + a)) /a/2; the 



: exp 



1 



V2xa) 



it2 



|0) , (35) 



where |0) is the vacuum state. If we set the par ame ter 



o2s 



2, i.e., the case of 3-dB squeezing, Eq. (34 1 is 



reduced to 



(x|sSCS±) = A4 



2 r 



±e 



(36) 
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That is, 



(x\ip n )=Af i _ ) nV- [e(*-^) + (_i)» e (»+v«) ] (37) 

with A/" ( _)n = 1/^/2 + 2(-l)»e- 2 ". 

Next, we can represent |0„) in the number-state basis 
with its wave-function defined in Ref. [8], 



/ 2 2n nl 
^{2n) 



(38) 



We use Eqs. (35 1 and (38) to obtain 




(39) 



where -ffn[x] is the nth-order Hermite polynomial. 



VIII. APPENDIX: B 



We can derive all the detailed expressions of the 
Wigner functions and the quantum fidelity in Sec. |III B| 
First, a Gaussian distribution of zero mean can be de- 
termined solely by its covariance matrix (CM) V = (Vij) 
where Vij — \{x-iXj + XjXi) — (xi)(xj) as 



W v (xi,...,x D ) = e 



-(l/2)x T ^- 1 x 



/ v / det(27rV0- (40) 



Here x = (xi, ...,xd) T is a D-dimensional vector, T de- 
notes the matrix transposition, and (x) is the expectation 
value of x. For example, a single-mode vacuum state can 
be represented by 



W yaD (fi) = Wv o (Pr,0i), 



(41) 



where its CM is given by Vq =1/4 with 1 being the 
2x2 identity matrix. Note that a coherent state with 
amplitude a is given by the same Wigner function with 
its center just displaced from the origin to a in phase 
space, i.e., W vac (/3 - a). 

We now describe the transformation of an input 
Wigner function to another under the operations such as 
single-mode squeezing, two-mode squeezing, and beam 
splitting. First, a single-mode squeezing transforms the 
phase-space variables as 



P^S^P where Si{s) 



e~ s 
e s 



(42) 



Here, s > (s < 0) indicates squeezing along the real 
(imaginary) axis in phase space. Accordingly, the covari- 
ance matrix V under the squeezing can be described by 
the transformation V — > Si(s)VSf (s) and the output 
Wigner function thus reads W Si ^ V()S t ^{P). 



Second, two-mode squeezing can be described similarly 
to the single-mode squeezing above by a transformation 
matrix 



5' 2 

with X = 



cosh s {I 

1" 

1 



1 X) +sinhs(A' 
and Z = 



(43) 



The Wigner function of a two-mode squeezed vacuum 
state in Eq. (22) can thereby be represented by 

W ad (f3 a ,f3 d ) = W Vad Wa,(3 d ) with 



V ad = S 2 {Vo®V )Sl 



(44) 



Third, a two-mode beamsplitting can be described by 
the transformation 



Pa 



-> B 2 {t) 



Pa 
Pb 



where -82(4) 



t -r 
r t 



Accordingly, two single-mode Gaussian states with their 
CMs, Va and Vb, are described under beam splitting as 

)(v a ®v b )bt (t) (Pa,Pb). (45) 

Tracing over one output mode after the beamsplitting 
of two Gaussian states gives a single-mode Gaussian state 
with a modified CM. Before presenting this formula, we 
will give some general formulas for integrating Gaussian 
functions. The first one is a partial integration of a Gaus- 
sian Wigner function that can be represented by 



d 2 P B W VAB {p A , p B ) = W Va(b) (p A ), (46) 



where Va<b) denotes the marginal CM for mode A, i.e., 
by removing block matrices related to mode B. The 
above result simply tells that the marginal distribution 
must be determined only by its original local CM. Note 
that the above and all the following results are not con- 
fined to the two-mode case but A and B modes can refer 
to any collective modes more than one. 

Next we present the partial integration of the product 
of two Gaussian functions 



d 2 (3 B W v (PA,PB)Wy(p B - a B ) 



(47) 



(v+o A t 

where O a denotes null matrix associated to mode A and 
V (V) refers to a CM of the two-mode (single-mode) 
state. Combining Eqs. (46) and (47), one can also easily 
derive another formula, 



d 2 p W v (p)Wy(p -a) = W v+V {a). (48) 



Now we are ready to obtain a Wigner-function formula 
associated with a beam splitter. Integrating Eq. (45) 
over variable P B using the formula (46), one can get 



W, 



d 2 p B W VA (tp A + rp B )Wv B (tp B - tPa) 

[B2(t)[V A <BV B )BT(t)] MB) {Pa) 



W t 2V A +r2V B {pA) 



(49) 
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If the above two Gaussian states have centers ua and 
ub other than the origin, the centers of the transformed 
Gaussian state becomes [B 2 (t)(aA, c<b) T ]a(b)- In the 
case of two multi-mode input Gaussian states, e.g., two 
two-mode states, the whole process of tracing out after 
beam splitting can be similarly represented by the trans- 
formation 

W VAB &v aB (Pa ,Pb,Pc,Pd)^ W Voat (Pa , Pb ) ■ (50) 
Here, the integration is performed over pc an d Pd and 

V mt = [B i {V A B®VcD)B{] AB(CD) 
= TV AB T T + RV cd R T , 
T = t{L © t 2 l, R = r\L © r 2 l, B A = 



T-R 
R T 



Note again that if there are nonzero centers 
(oca, &b, «Ci chd) T in the original Wigner function, 
those transform to [B±(aA, as, etc, &d) T ]ab(cd)- 

Using the general results for Gaussian integration so 
far, analytical expressions for all the Wigner functions 
in the main text are readily obtained. First, the Wigner 
functions in Eqs. (23) and (24) read 



W a (p a ) 



1 

Pd 



R, 



-It I d 2 p d Wv ad Wa,f3d)W Vo (/3 d ) 
1 

1 — sech 2 s = tanh 2 s. 



[W cosh2sVo (p a ) - sech 2 sW Vo (p a )] , (51) 



(52) 



Next, for the state in Eq. (|25|), we obtain 

W ac (P a ,Pc) 



1 

Pd 



W S2 l( cosh2s V )(BVo]ST Wa, Pc) 



sech 2 sW S2[VoS)Vo]S T(f3 a ,f3 c ) 



1 

Ph. 



'W v . e (Pn,Pc) ~ sech 2 sW Vad (p a ,Pc)], 



with V ac — ^[(cosh 2 sVq) © Vo]S 2 . Beam splitting these 
states gives the Wigner functions in Eqs. (26) and (27), 



W abc (Pa,Pb,Pc 



Pd 



[W vm (P a ,Pb,Pc) 



-sech 2 sW" v ( 2 ) (P a , p b , p c )] , 



W' abc (P a ,P b , p c ) = — [W V , W (P a , p b , P c ) 



p, 



-sech 2 sW" y / ( 2) (P a , p b , p c )] , 



where 



y^c^Bab[Vac®(V )b\Blb 



K (2) = B nh V„,B 



b^adOab, 



,(2) D T 
bci 



T/'W _ D T/ 1 rT t/(2J _ r> 1/ 2 R 
V abc - B bcV a bc U bci V abc - B bcV ab c U i 

B ab = B 2 (t 1 )®X, and B bc =X®B 2 (t). 



Finally, on-off detecting on modes b and c gives the out- 
put Wigner functions as 



W out (p a ) = — ^- J2 \A^W v aMPa) 

i^d^bc , L ° ut 

-sech 2 s A^W vl 2 :3) (pa) 
4 

Pbc 



J ') - sech 2 s 



3=1 



where 2 = 1,2 and 



out. 
(<,2) 



V 



/(i) 
a(6c) 



[^g)c + Oa©(^)c]^ 



T/(<.3) _ 



(M) 



^a6c + °* © (%)& ® (Vo)c]&, 



A(^) = l, 
A (^,2) = _ 

^(i,3) = _ 

A (i,4) 



Wdet[V$> (c) +V ] 



Here, the symbol C denotes Schur complement [BU] with 
the notation 

(Vab)a = Vb,b - Vb,aVa^ a Va,b, 

(Vab) C b = V a ,a - V a ,bV b ] b V b ,a, (53) 

where Va.a, Va,b and so on are block matrices in 
the whole matrix Vab associated to the corresponding 
modes. 



The fidelity in Eq. ( 30 ) can also be calculated straight- 



forwardly as the Wigner function of the target sSCS in 
Eqs. (31 ) and (32) can be represented by Gaussian func- 



tions, i.e., 

W^scs 63) = M 2 [W Vo (P r - a, Pi) + Wvo (P r + a,Pi) 

-2W Vo (P)cQs(4api)], (54) 

W sSCS (P) = [W scs (P)]^X, )VoST(s , y (55) 



Using these expressions and the formula in Eq. ( 48 ) , we 
obtain the fidelity 

-sech 2 S A^W [v ^ ])+Si{s , )VoST(sl)] (0,iae- s ') 
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